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Abstract. We study the quasineutral limit of the isothermal Euler-Poisson system describing a 
plasma made of ions and massless electrons. The analysis is achieved in a domain of R'^ and thus 
extends former results by Cordier and Grenier [Comm. Partial Differential Equations, 25 (2000), 
pp. 1099-1113], who dealt with the same problem in a one-dimensional domain without boundary. 

1. Introduction 

We consider the isothermal Euler-Poisson system which models the behaviour of ions in a back- 
ground of massless electrons. The massless assumption implies that the electrons follow the classical 
Maxwell-Boltzmann relation: denoting by Ue their density and after some suitable normalization 
of the constants, this reads 

rie = e~\ 

where —cj) is the electric potential (auto-induced by the charged particles), which satisfies some 
Poisson equation. 

Introducing the characteristic observation length L and the Debye length Xd = yf^^^^ the 

plasma (where T* is the average temperature of the ions, and Ni their average density), we are 
interested in the study of the behaviour of the system, when 

which corresponds to the (usual) physical situation where the observation length is much larger 
than the Debye length (indeed, in usual regimes. Ad ~ 10~^ — lO^^m). Loosely speaking, the 
Poisson equation satisfied by —cj) then reads: 

e'^Ax(j) = rii - Ue- 

In situations where e <C 1, this indicates that the plasma can be considered as being "almost" 
neutral (that is iii ~ ng): hence the name quasineutral limit for the limit e — ?■ 0. However, in a 
domain with boundaries, it is well-known that there can be some interaction between the plasma 
and the boundary, so that small scale effects are amplified and boundary layers generally appear. 
This implies that quasineutrality breaks down near the borders. 

In this paper, our goal is to rigorously study this phenomenon of boundary layers by investigating 
their existence and their stability. More precisely, the system we work on is the following dimen- 
sionless isothermal Euler-Poisson system, for t > and x = [xi,X2,x^) = {y, x-^) G M.^ := M? x M+: 

' dtn + div (nn) = 0, 

(1.1) } dtu + u-Vu + T V ln(n) = V(/>, 

, e^Acf) + e-* = n, 

where n is the density of ions, u = {ui,U2,u^) = {uy,u^) is their velocity field, and —cj) is the 
electric potential. 
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We complete this system with the bomidary conditions 
(1-2) U3\x3=0 = 0, 01x3=0 = (l^b, 

where (pb is some prescribed potential. We assume that (pf, = (j)b{y) is smooth and compactly 
supported. Otherwise, it is also possible to consider with minor modifications the case where (pb is 
independent of y, which corresponds to the perfect conductor assumption: 

Vy(t> = 0. 

The boundary condition on u corresponds to a non-penetration condition. From the physical point 
of view, this can be interpreted as a kind of confinement condition. Non homogeneous boundary 
conditions are also physically pertinent, and some cases could be treated in our framework: we refer 
to Remark [3] and Section 15. 2i Also we shall focus in our analysis on the three-dimensional case 
and the simplest domain Mi]_. Nevertheless, our analysis can be easily extended to any dimension 
or more general domains. 

The Euler-Poisson system and its asymptotic limits have attracted a lot of attention over the 
past two decades, possibly due to its many applications to plasmas and semiconductors (see for 
instance [26^ I28j). Among many possibilities, the so-called quasineutral limit, zero-electron-mass 
limit, zero-relaxation-time limit have been particularly studied. 

For the sake of brevity, we only focus here on the quasineutral limit and describe some math- 
ematical contributions. For the corresponding Euler-Poisson system describing electrons (with 
fixed ions), in a domain without boundary, there have been many mathematical studies of the 
quasineutral limit (for various pressure laws) and the situation is rather well understood; we refer 
for instance to [24:\ [321 [271 l^Ol [33] and references therein. The quasineutral limit for Two-fluid 
Euler-Poisson systems (one for ions and one for electrons), still without boundary, has been also 
recently investigated in [23l [22] . 

Considering the dynamics of ions with electrons following a Maxwell-Boltzmann law (which is 
the one we are particularly interested in), Cordier and Grenier studied in [6] the quasineutral limit 
of the Euler-Poisson system (|l.ip . in the whole space M. A kinetic version of the problem was 
also recently studied in [2^. We also mention the work [19| of Guo and Pausader who recently 
constructed global smooth irrotational solutions to (|l.ip with e = 1 in the whole space M^. 

For quasineutral limits in presence of boundaries, there have been very few rigorous studies, at 
least to the best of our knowledge. When the Poisson equation is considered alone (the density 
of the ions being prescribed), a boundary layer analysis was led by Brezis, Golse and Sentis [1] 
and Ambroso, Mehats and Raviart [2\. For the Poisson equation coupled with some equations 
for the plasma flow, the are only two main results we are aware of. The first one is due to 
Slemrod and Sternberg [36], who dealt with the quasineutral limit of some Euler-Poisson models 
(with massless electrons) in a one-dimensional and stationary case (see also Peng [31j for general 
boundary conditions). The second is due to Peng and Wang [30] and Violet [39], who studied the 
quasineutral limit of some Euler-Poisson model (with fixed ions) for stationary and irrotational 
flows. 

To conclude this short review, let us finally mention some other interesting directions which are 
related to the issue of the quasineutral limit in a domain. In [37], Suzuki studies the asymptotic 
stability of specific stationary solutions to some Euler-Poisson system, which are interpreted as 
"boundary layers" (see also the paper of Ambroso [1] for a numerical study). In [201 [TO] and 
references therein, Slemrod et al. study the problem of formation and dynamics of the so-called 
plasma sheath. We refer also to the work [9] of Crispel, Degond and Vignal for a study of a two-fluid 
quasineutral plasma, based on a formal expansion. 
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(2.1) 



In this paper, we study the quasineutral hmit of ()l.ip - (jl.2p . As far as we know, our results 
which are described in the next section, give the first description and stabihty analysis of boundary 
layers for the quasineutral limit for the full model in a non stationary setting (and we are in a 
multi-dimensional framework as well). 

2. Main results and strategy of the proof 

By putting formally e = in the Poisson equation of , we first directly obtain the neutrality 
condition: 

n = e-'^, 

so that one gets from (jl.ip in the limit e — )• the following isothermal Euler system: 

{dtn + div (nu) = 0, 
dtu + u-Vu + (T* + l)Vln(n) = 0. 

It is natural to expect the convergence to this system. The previous system is well-posed with the 
only boundary condition 

U3\j:^=0 = 0. 

Because of the boundary conditions (jl.2p . we would like to satisfy also the condition 

(p\x3=0 = 4>b- 

Nevertheless, the solution (n, n) of (j2.ip cannot in general satisfy in addition n\xg=Q = e~'^^ and 
hence, we expect the formation of a boundary layer in order to correct this boundary condition. 

We shall consider solutions of ()2.ip with no vacuum: we fix a reference smooth function rf^^ 
which is bounded from below by a positive constant and consider solutions of ()2.ip for which n 
is bounded from below by a positive constant and such that n — rf^^ is in . The main result 
contained in this paper is a rigorous proof of the convergence to the isothermal Euler system (j2.ip . 

Theorem 1. Let {n^,v?) a solution to (HH) such that (n^ - n'^^'f ,u^) G C°([0, T], i?''(M3 )) with 
s large enough. There is > Q such that for any e £ (0,eo]; there exists [n^ ,u'^) a solution to 
()l.ip - ()1.2p also defined on [0, T] such that 

(2.2) sup (^||n^ - n°||^2(R3 ) + \\u^ - u°||l2(r3_)) ^e-^o 0. 

Furthermore, the rate of convergence is 0{^/e). 

Remark 1. Let us recall that following the classical theory of the initial boundary value problem 
for the compressible Euler equation (see \35\ I18j for example), asking {nP,u^) to be a smooth so- 
lution to the compressible Euler system imposes some compatibility conditions on the initial data 
(n°|t=o,ii°|t=o)- 

Boundary layers are hidden in the statement of Theorem [1] (recall that they are expected to 
appear since the boundary conditions do not match). Nevertheless, in order to give a proof, we 
will have to describe them precisely. Indeed, for the system (jl.ip in multi-d, the only existence 
result that is known is the local existence result of smooth (i/™, m > 5/2) solution that can be 
obtained from the theory of the compressible Euler equation (see j35l [18] for example) . A nontrivial 
part of the statement of Theorem [1] is thus that the solution of p.ip is defined on an interval of 
time independent of e. In this is a consequence of the fact that (jl.ip enjoys estimates that 
are uniform in e. These estimates were obtained in [6] by using that (jl.ip can be rewritten as a 
symmetric hyperbolic system perturbed by a large skew-symmetric operator and thus can be cast 
in the general framework of [13] that generalizes the classical theory of Klainerman and Majda 
[25] on the incompressible limit. When boundaries are present it seems unlikely to derive directly 
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uniform estimates of (jl.ip . Indeed the expected description of the solution {rf,u^,(f)^) of (jl.ip 

is under the form 

(2.3) 

= n°(t, x) + iV° (t, y, ^) + 0(e), = </)°(t, x) + $° (t, y, ^) + 0{e), = u^{t, x) + 0{e) 

where {n^,u^) is the solution of (|2.ip . = — logn^ and the profiles A^^, (the boundary layer 
profiles) are smooth and fastly decreasing in the last variable and are added in order to match the 
boundary condition ()1.2p . In particular, we require 

cf>^{t,y,0) + ^^ {t,y,0)=4>'{y). 

Note that at leading order, there is no boundary layer for the velocity but that in general, since 
the solution of ()2.ip does not satisfy n\x^=o = e~'^'', the profiles 4>*^ and are not zero for positive 
times. This yields that cannot satisfy uniform if™, m > 5/2 estimates. 

In order to prove Theorem [H we shall use a two step argument. We shall first prove that, up to 
any order, there exists an approximate solution (n^^^, ti^pp, 4>app) for in the form of a two-scale 

asymptotic expansion, and then, we shall prove that we can get a true solution of (ll.ip defined 
on [0, T] by adding a small corrector to the approximate solution. The description of the solution 
that we get is thus more precise than stated in Theorem [TJ We refer to Theorem [3] and Corollary 
[T] which are stated in Section [H In particular, as a simple corollary of our estimates, we will get 
that the behavior in is given by: 



+ \\U^ — u'^ll rooms -1 



sup 

(2.4) ^ 



'■-) 



0. 



This two step approach has been very useful to deal with other asymptotic problems involving 
as above boundary layers of size e and amplitude 0(1) like the vanishing viscosity limit to general 
hyperbolic systems when the boundary is non-characteristic [12l [151 IIZl EH] > see also [38] , or the 
highly rotating fluids limit [16[ \TT\ I34j . A common feature of these works is the fact that the 
true solution of the system remains close to the approximate one on an interval of time that does 
not shrink when e goes to zero (i.e. the stability of the approximate solution) requires a spectral 
assumption on the main boundary layer proflle (that is automatically satisfled when the amplitude 
of the boundary layer is sufficiently weak). One quite striking fact about the quasineutral limit 
is that we shall establish here the existence and unconditional stability of sufficiently accurate 
approximate solutions involving boundary layers of arbitrary amplitude. 

Since the boundary layer profiles solve ordinary differential equations, we shall be able to get 
their existence (without any restriction on their amplitude) by using classical properties of planar 
Hamiltonian systems. In order to prove the stability of the approximate solution, there are two 
difficulties: the fact that the electric field is very large when e is small (this difficulty is already 
present when there is no boundary) and the fact that the linearization of (jl.ip about the approxi- 
mate solution contains zero order terms that involve the gradient of the approximate solution and 
thus that are also singular when e goes to zero because of the presence of boundary layers. We shall 
avoid the pseudo-differential approach of [6] which does not seem to be easily adaptable when there 
is a boundary and get an type estimate which is uniform in e by using a suitably modulated 
linearized version of the physical energy of the Euler-Poisson system. We refer for example to [14[ [5] 
for other uses of this approach in singular perturbation problems. In order to estimate higher order 
derivatives, we shall proceed in a classical way by estimating conormal derivatives of the solution 
and then by using the equation to recover estimates for normal derivatives. Since the boundary 
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is characteristic this does not provide any information on the normal derivative of the tangential 
velocity but we can estimate it by using the vorticity equation. 

Theorem [1] (as well as Theorem [3] and Corollary [1]) can be generalized to more general settings; 
this is the purpose of the two following remarks. 

Remark 2 (General domains). Less stringent settings for the domain could also be treated. Instead 
of the half-space M5|_, we can as well consider the case of a smooth domain $7 C M'^ whose boundary 
is an orientahle and compact manifold. We refer to Section \5.1\ for some elements on how one 
could adapt the proof. 

Remark 3 (Non- homogeneous boundary condition on u). Finally, it is possible to adapt the proof 
in order to study, instead of the non-penetration condition on u, the case of some subsonic outflow 
condition on the normal velocity u^. We refer to Section \5.^ for some details on what changes in 
the proof. 

This paper is organized as follows. Section [3] is dedicated to the construction of a high order 
boundary layer approximation to a solution of (jl.ip (see Theorem [2]) . In Section U we prove the 
stability of this approximation ( c/ Theorem [3]) , thus proving the convergence to isothermal Euler. 
The possible extensions and variations about the result are discussed in the last section. 

3. Construction of boundary layer approximations 
We look for approximate solutions of system (II. ip of the form 

K 

ina,Ua,(t>a) = ^ {u' {t , x) , u\t , x) , {t , x)) 

.E.'(iv(v.f).^-(*,.?),.'0,.?)) 

i=0 

where K is an arbitrarily large integer. The coefficients (n*,u*,(/)*) of the first sum depend on x. 
They should model the macroscopic behaviour of the solutions. The coefficients {N^, U^, $*) of the 
second sum depend on t,y, but also on a stretched variable z = ^ G M+. They should model a 
boundary layer of size £ near the boundary. Accordingly, we shall ensure that 

(3.2) {N\U\<^')^0, as z^+oo 

fast enough. In order for the whole approximation to satisfy (jl.2p . we shall further impose 
u\{t, y, 0) + Ui{t, y, 0) =0 for all i > 0, 

/(t,y,0) + $°(t,2/,0) = 0^(t,y,O) + ^'{t,y,0) = for alH > 1. 

Thus, far away from the boundary, the term {n^,u^,(l)^) will govern the dynamics of these approx- 
imate solutions. As explained in the previous section, we expect {n^,u^) to satisfy the following 
compressible Euler system 

{dtu + div (nu) = 0, 
dtu + u-Vu + {T + l)Vln(n) = 0, 



(3.4) 



together with the relation = e . As we will see below, it will also satisfy the non-penetration 
condition 

(3.5) U3\x3=o = 0. 



Our results are gathered in the: 
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Theorem 2. Let K G N*,m G N* such that m > 3. For any data (riQ,Ug) satisfying some 
compatibility conditions on {xa = 0} and such that (ti-q — n''''^-^ , Uq) G H"^'^'^^~^^{M.^), there is T > 
and a smooth approximate solution of (jl.ip of order K under the form p.ip such that 

i) (nO-n''^-^,'u°) G C°([0,T],F'"+2^+3(m3 j)^ ^0 = _iognO, (n^u") is a sotoion to the isother- 
mal Euler system (|3.4|) wii/i initial data {n^^u^). 

n)Ml<i<K, {n\u\(t)') G C°([0, T], i/™+3(M^)). 

iii) VO < i < K, N^,U^,^^ are smooth functions that belong together with their derivatives to 
the set of uniformly exponentially decreasing functions with respect to the last variable. 

iv) Let us consider (n^, u^, (/>^) a solution to and define: 

n = -Ua, U = U'^ -Ua, 4> = - 4>a- 

Then (n, n, (p) satisfies the system of equations: 



(3.6) 



' dtn + {ua + u) ■ Vn + n div (n + Ua) + div (uau) = e Rn, 

\na + n Ua \na + n' 



^ e^Acj) = n- e-<^° {e-* - l) + e^+^R^. 
where Rn, Ru, R(i> cti"^ remainders satisfying: 

(3.7) sup \\V^Rn,u,cj,\\L2(R^) < C'ae~°^ Va = (ai,Q;2,a3) S N^, \a\ < m, 
[o,T] + 

with Ca > independent of e. 

The rest of this section is devoted to the proof of this theorem. 

3.1. Collection of inner and outer equations. First, we derive formally the equations that 
the leading terms of the expansion (j3.ip should satisfy, to make it an approximate solution of 
(jl.ip - (|1.2p . The solvability of these equations will be discussed in the next subsection. 

The starting point of the formal derivation is to plug the expansion (|3.ip into (jl.ip . By consid- 
ering terms that have the same amplitude (that is the same power of e in front of it), one obtains 
a collection of equations on the coefficients {n^,u^,(j)^) and (A^*, [/*, $*). As usual, we distinguish 
between two zones: 

(1) An outer zone (far from the boundary), in which the boundary layer correctors are ne- 
glectible. 

(2) An inner zone (in the boundary layer), in which we use the Taylor expansion 

n^{t,y,x^) = n^{t,y,ez) = n^{t,y,0) + edsn'''{t,y,0)z + ... (the same for n* and 0*) 

to obtain boundary layer equations in variables {t,y,z). For clarity of exposure, for any 
function / = /(t, x), we will denote by Tf the function (t, y, Z) i— f{t, y, 0). 

In the outer zone, by collecting terms that have 0(1) amplitude, one sees that [n^,u^) satisfies 
([TT]) . and that n° = e"*^". 

In the inner zone, by collecting terms that have 0{e~^) amplitude, one has the relations 

r d4{Tn'> + N'){Tu'i + Ui)) = 

\ (Tul + C/3°) d,Ul^ + T'd, ln(rn° + N^) = 9^$° 

The first line implies that (Fn^ + A^'')(rti3 + f/g) does not depend on z. From condition ()3.3p . we 
deduce 

(3.8) (rn° + (Tul + Ul) = 0. 



Now, the second line also reads 

(3.9) (rn° + iv°) {Tul, + u^) d.ul, + ra,iv° = (rn° + iv°) 9,$°. 

Combined with the previous equation, this gives 

Td^N^ = (rn° + N^)d,^^ 

or equivalently (still using (|3.2p ) 

(3.10) (iV" + rnO) = rn°e*°/^' 

Back to dSl), we obtain that Tul + ^^3 = 0' and by ([221), both = and rn[| = 0. We recover 
as expected that {n^^u^) satisfies formally (j3.4p . together with the non-penetration condition ([33]), 
and the relation vP = e~^'' . 

Still in the inner zone, collecting the 0(1) terms in the third line of (jl.ip . we obtain 
As nP = e~'^° , this last equation can be simplified into 
Together with (|3.10p . we end up with 

(3.11) a2«>o + rn°5(«>0) = 0, S{<^>) := e-* - e*/^\ 

Note that this is a second order equation in z, so that with the boundary conditions 

(3.12) $0|^=o = </'°U3=o, ^'°U=+oo = 

it should determine completely Indeed, the existence of a unique smooth solution of (|3.1ip - 
(j3.12p with rapid decay at infinity will be established in the next subsection. As is determined, 
so is iV° by relation (l3J0]l . 

This concludes the formal derivation of {nP , , (j)^) , and {Ul^,^^,N^). We will now turn to the 
equations satisfied by {n^ and (C/3 , Uy,^^, N^). More generally, we will derive for alH > 1: 

(1) a set of equations on (n*, u*, </>*), whose source terms depend on the coefficients {rJ^ ,<f)^) 
and their derivatives, k < i — 1. 

(2) a set of equations on ([/I, C/*~^, A^*), whose sources depend on (the trace at 2:3 = 
of) {n^ and their derivatives, /c < i, as well as on [U^ ,Uy~^ , N^) and their 
derivatives, k < i — 1. 

Indeed, in the outer zone, collecting terms with amplitude 0(e*) in (jl.ip leads to systems of the 
type: 

' dtn' + div (n%^) + div (n*n°) = 



(3.13) 



^0 



, - cP' = n' + 

where the source terms f^ucf) depend on {n^,u'',(j)^) and their derivatives, k < i — 1. 

In the inner zone, collecting 0(e*) terms in the Poisson equation of (jl.ip . we get: 
(3.14) ^2$^ - rn^e"*" ($* + T(j/) = {N' + rn') + Fl 



where F|, depend on {n'',(f)^) and on {N'',^^), k < i — 1. Then, collecting the 0(e* ^) in the 
equation for ion density n, taking into account that = and Tu'^ = 0, we get 

(3.15) a, {{N^ + rn°)(C/^ + rn^)) + div, ((iV" + Tn'){W-' + T^-^)) = F^. 

The equation on Uy yields 

dtiU'y + r^°) + iU^ + + Tds4) d^U^y 

+ + ru°) • (C/° + rn°) + (r + l)Vy In(rnO) = 

and for i > 2: 

Finally, the equation on us yields 

(3-18) ro. [wTf^) = + ^3- 

Again, the source terms F^, Fy and Fg depend on {n^ ,<j)^), and on {U-^ ,Uy~^ , N^), for 
indices k < i — 1. Note that 1)3.20 requires the following compatibility conditions: 

— n°(/)*|^3=o = '^*U3=o + -^4U=oo, 
divy (n°ny"^) |a;3=o = -FAfU=oo, 

-^3^=00 = 0. 



(3.19) 



It remains to show how to deduce the profiles from the equations. First of all, we can derive Uy 
from (I3.15P (with i = 1) and (I3.16p . We integrate (I3.15P from to z, and with (13. 3|) deduce 

(3.20) {Ul + Tui) = - ^olrn^ [ i^^^v ((^° + ^^'Wv~' + ^4'')) " ^n) 

For 1 = 1, this last relation allows to express + Tu\ in terms of Uy + Fii^ (and coefficients of 
lower order). One can then substitute into (|3.16p to have a closed equation on Uy. We notice that 
it has the trivial solution Uy = 0. In all what follows, we shall restrict to approximate solutions 
([3T]) satisfying f7° = 0. 

As [/^ = 0, the equation ()3.17p simplifies into 

(3.21) dtUi-^ + Fn° • VyU'y-^ + U^^^ ■ VyVul = W 
with 

W := Fi - {dtTv^y^ + FnO • VyVv"-^ + Tv'-^ ■ Vj,Fn°) . 

One can now derive recursively (n*, u*, </»*) and (C/|, f/*^^, A'^*), i > 1. We take i > 1, and 
assume that all lower order profiles {n^ , , cj)'') and (C/3 , A"^), k < i — 1 are known. 

From the equations satisfied by the {n'^,u^,<p^ys, the last three compatibility conditions in (j3.19p 
are satisfied. In particular, the function in (|3.2ip decays to zero at infinity. Thus, the linear 
transport equation (I3.2ip yields a decaying solution Uy~^ (with the special case C/° = when 
i = 1). Then, one obtains from (j3.20p the expression of C/3 + rii|. From the decay condition p.2p . 
we deduce that 

(3.22) ^1^3=0 = / (div, ((AtO + FnO)(C/-i + F^.-!)) - F^) 

n 1x3=0 JO 



This boundary condition goes with the hyperbohc system on (n^,u^) (see (j3.13p ). Together with an 
initial condition which is compatible with the boundary condition, it allows to determine (n*, u*, (/>*). 
This will be explained rigorously in the next subsection. Note that, as (n*, u*, solves p.lSp . the 
remaining compatibility condition in (|3.19p is satisfied. Using again (|3.20p . one eventually gets U^. 

To end up this formal derivation, it remains to handle and A^*. Combining (13.14p and (|3.18p 
(together with the decay condition (|3.2p ) leads to a second order equation on 

(3.23) af^* + rn°5'($°) = G$ 

where we remind that S{^) = e~* — e^^'^' and is a source term that depends on the lower order 
terms, decaying to zero as z goes to infinity. Again, the solvability of this system with Dirichlet 
conditions 

(3.24) ^%=o = -<^ix3=o, $i.=+oo = 

will be established later on. Once <1>* is known, is determined through (j3.18p . 

3.2. Well-posedness of the reduced models. To complete the construction of the approximate 
solutions and get Theorem [21 we must establish the well-posedness of the inner and outer systems 
derived in the previous paragraph. 

The well-posedness of the inner systems is classical. Indeed, the system (j3.4p - (|3.5p on {nPjU^) 
is a standard compressible Euler system. For well-posedness, we consider an initial data of the 
form (ng = e'^o^u'^), with («[], ng - n''^-^) E i/'"+3+2^(M3 )! Following [iSlllIj, if this initial data 
satisfies standard compatibility conditions at the boundary {x^ = 0}, there exists a unique solution 
(nO = e-'^°,'u°), with 

(nO - n''^^,M°) E C7°° {[0,T]; H'^+^+'^^ {RD^) , for some T > 0. 

As regards the next systems (|3.13p - (l3.22l) . they resume to linear hyperbolic systems on (n*,n*), 
z > 1, with a source made of some derivatives of previously constructed {n^ , (j)^), j < i. For 
initial data (tiqjUq) E H^^^'^'^^ ""^^ (M.^)'^ , chosen in order to match compatibility conditions they 
have again unique solutions 

{n\u') E C°° ([0,T];i?'"+3+2X-2i(^3 ^ 

The "loss" of regularity is due to the derivatives in the source. Likewise, one shows that: 

W E C~ ([0,T];if™+3+2^-2*^]^3 ^ 

The main point is the resolution of the nonlinear boundary layer system p.lip - (j3.12p . Note 
that t and y are only parameters in such a system, involved through the coefficient 'ynP and the 
boundary data. For each {t,y), ()3.1ip is an ordinary differential equation in z, of the type 

(3.25) $" + 7cS($) = 0, 7>0, 

and we want to prove that, for any constant 0, it has a (unique) solution $ that connects (j) to 0, 
with exponential decay of $ and its derivatives as z goes to infinity. 
Therefore, we rewrite (j3.25p as a Hamiltonian system, 

(3.26) ^ Q = V^Hip,-^), p := Hip,<^) := y + r($) 
where T is an antiderivative for 75. We choose the one that vanishes at 0: 

r(cl>) := -7 (e~* + re*/^') + 7(1 + T'). 
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The Hamiltonian H is of course constant along any trajectory. The Unearization of (j3.26p at the 
critical point (p, = (0, 0) yields 

Hence, we get that (0, 0) is a saddle fixed point, and using the stable manifold theorem that its 
stable manifold is locally a curve, tangent to ^ Vii^+^T^) ^ . Moreover, as H{0,0) = 0, the stable 
manifold is exactly the branch of {H{p, <I>) = 0} given by the equations 



p = _^y^e-* + T^e*/"^' - 1 - for $ > 0, p = /rVe"* + T^e*/"^' - 1 - for ^> < 0. 

In addition, any solution starting on this branch decays exponentially to 0, and by using the 
equation (I3.25p . so do all its derivatives. Finally, one needs to check that for any tp, this branch 
has a unique point with ordinate V- It is indeed the this branch gives p as a decreasing 

function of going to infinity at infinity. This concludes the resolution of the o.d.e, and yields in 
turn the solution <I>^ of our boundary layer system. Its regularity with respect to (t, y) follows from 
regular dependence of the solution of (13.250 with respect to the data. 

Last step is the well-posedness of the boundary layer systems (|3.23p - (|3.24p which allow to de- 
fine <I>* and A^*. Again, t,y are parameters, and the point is to find a (unique) solution to the 
inhomogeneous linear ODE 

(3.27) ^" + 75'(¥)«> = F, 

with exponential decay at infinity and a Dirichlet condition $|^=o = V'- Here, $ and F are arbitrary 
smooth functions of z, decaying exponentially to 0. 

Note that, up to consider ^{z) = ^{z) — ipxi^)^ with x £ ^"^(1^+) satisfying x(0) = 0, we can 
always assume that 1^ = 0. After this simplification, observing that S"(<I>) < 0, we can apply the 
Lax-Milgram Lemma: it provides a unique solution ^> G F(J(1R+) of ([327]) • By Sobolev embedding, 
it decays to zero at infinity, and it is smooth by elliptic regularity. It remains to obtain the 
exponential decay of ^, hence we write (I3.27P as 

^" + 75'(0)$ = G, G:= {-fS'{0) - -fS'^)^ + F 

Using the Duhamel formula for this constant coefficient ode, and the exponential decay of G, one 
obtains easily that any bounded solution decays exponentially. 
Finally one deduces recursively that W also decays exponentially. 

Summing up what has been done in this section, we have proved Theorem [2j 

4. Stability estimates 

This section is devoted to the stability of the boundary layer approximations built in the previous 
section. 

Let us write the solution (n^,n'^,0^) of (jl.ip under the form 

= Ua + n, = Ua +U, (f)'^ = (t)a + 4> 

where na,Ua,4>a are shorthands for n'^ppjU^ppjCpapp (defined in (j3.ip ). Then, we get for {n,u,(j)) the 
system 



(4.1) 



' dtn + {ua + u) ■ Vn + n div {u + Ua) + div {uau) = e Rn, 

\na + n Ua \na + n; 
, e'^AcP = n- e-^'^ (e""^ - l) + e^+^i?^. 
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together with the boundary conditions 
(4.2) ^31x3=0 
and the initial condition 
(4.3) 



1X3=0 



u\t=o 



e-^+^uo, n\t=o 



Observe here that in (j4.ip . Rn, £^ Ru and e^^^Rfj) are remainders that appear because (n^, Ua, (pa) 
is not an exact solution of (jl.ip . 



The main result of this section is 

Theorem 3. Let m > 3, and {jiq^uq) E i/™(R5|_) some initial data for (j4.3p . satisfying some 
suitable compatibility conditions. Let K S N*,if > m and {na,Ua,4>a) 0,^ approximate solution 
at order K given by Theorem [H which is defined on [0, Tq] . There exists £q such that for every 
e G (0,60]; ihe solution of (|4.ip - (|4.2p - (j4.3p is defined on [0,To] and satisfies the estimate 

el"l||0"(n,n,(/>,eV(/>)||i2(K3) < Ce^, Vt G [0,r°], Va G N^ |a| < m. 

One can then remark that as a simple rephrase of Theorem [31 we obtain: 

Corollary 1. Let m > 3. Let K G fi*,K > m and {na,Ua,4>a) an approximate solution at order 
K given by Theorem\^ which is defined on [0,To]. There exists such that for every e G (0,eo]; 
there is a solution [n^ ,u'^ ,(f)^) to (II. ip which is defined on [0,To] and satisfies the estimate: 



(4.4) 



n 



Ua.U 



< Ce 



G [0, r° 



Ln particular, we get the L^ and L°° convergences as e ^ Q: 



(4.5) 



sup 

[0,To] 

sup ( 

[0,To] ^ 

sup 

[O.To] 



n 



1%) 



+ \\n 



0, 



+ M 



+ 



0. 



Of course, this contains Theorem [H 



Prom now on, our goal is to prove Theorem [3l 

For e > fixed, since in the equation (j4.ip . the term involving cf) in the second equation can 
be considered as a semi-linear term, the known local existence results for the compressible Euler 
equation ([35l[18] for example) can be applied to the system (j4.ip . Let us assume that (no,uo) S 
H"^(R.^) for m > 3 and that it satisfies suitable compatibility conditions on the boundary, then 
there exists T^ > and a unique solution of ()4.ip defined on [0,T'^] such that u G C{[0,T^), H"^) 
and that there exists M > satisfying 
(4.6) 

na{t,x)+n{t,x) > l/M, e""^" (l+min(/io, ^1)) > \\x{ua+u)3\\L-- < ^3/42^, Vt G [0,T^) 

where 

e--^ - 1 + c' 



(4.7) 



/io(0) := 



hi{<P) :-- 



1 
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and xi^s) •= xi^s/^) where x is a smooth compactly supported function, equal to 1 in the vicinity 
of zero and 5 > is chosen so that 



|x(^3M)((na + n)3)|t=o||L- < 7^72. 



Note that we can always chose 6 in this way since at t = we have [{ua + u)^^ |j=o,x3=o = 0. 

The difficulty is thus to prove that the solution actually exists on an interval of time independent 
of e. We shall get this result by proving uniform energy estimates combined with the previous local 
existence result when the initial data and the source term are sufficiently small (i.e. when the 
approximate solution {na,Ua) is sufficiently accurate). Note that there are two difficulties in order 
to get useful energy estimates. The first one is the singular perturbation coming from the Poisson 
part, the electric field cannot be considered as a lower order term uniformly in e. The second one 
comes from the boundary layer terms contained in the approximate solution which create singular 
terms, for example, in the second line of ()4.ip . the zero order term ^ is singular in the sense 
that when (j4.6p is matched, we only have the estimate 



Vua n 



Ha ria + n 



L2 



< lll II 



The main step will be the obtention of uniform estimates for "linearized" systems. 

4.1. Energy for the quasineutral Euler- Poisson system without source. We start by re- 
calling that the isothermal Euler-Poisson system without source (which corresponds here to the 
case (ph = 0) has a conserved physical energy, which is given in the 

Proposition 1. Let e > and {n,u,4>) a strong solution to (jl.ip on [0,T] with (pb = 0. We define 
the energy functional: 

(4.8) Seit) ■■= I [ n\u\'^dx + T [ n{logn-l)dx+ [ (1 - 0)e-'^dx + ^ / \V^^\'^dx. 

^ Jul JRl JrI ^ Jul 

Then for any t £ [0,T], S^it) = <?e(0). 

This property will never be used in the sequel; nevertheless, the and higher order stability 
estimates which follow are obtained via the study a modulated version of this energy. Thus, for 
clarity of exposure, we briefiy present the proof showing that the energy is conserved, as it is much 
easier to follow but share the same spirit with the subsequent ones. 

Proof. We compute the derivative in time of the first term of £, by using the transport equation 
satisfied by n in (11. ip (which corresponds to the convervation of charge): 

— [ n\u\'^dx = [ dtn\u\'^dx + 2 f nu ■ dtudx 

at Jr3_ 7r3^ Jj^3^ 

= — V ■ {nu)\u\'^ dx + 2 / nu-dtudx. 
JrI JrI 

Using the non-penetration boundary condition for u, we have by integration by parts: 

V ■ {nu)\u\'^ dx = / (nu) ■ V\u\'^dx = / {nu)u-Vudx. 



Now thanks to the equation satisfied by u in (jl.ip . we can write: 

/ nu ■ dtudx = — / {nu)u ■ Vudx — / nn-Vlogn+ / nu-Vcpdx. 

Jr^. Jr3 
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By the equation satisfied by n and the non-penetration condition, we infer, by another integration 
by parts, that 

T* / nn-Vlogn=-r' [ dtnlogn = - —T^ [ n{logn - l)dx. 

Likewise, we get, using once more the conservation of charge, 

nu • V(j)dx = — V • {nu)(j)dx 



I 



dtn (pdx. 

Deriving with respect to time the Poisson equation in we obtain 

nu ■ V4>dx = / e"^ dtA(j)(j)dx + / dte~'^(j)dx. 



The last term of the r.h.s. is treated exactly like the pressure term. Considering the first one, by 
integration by parts, we obtain: 

dtA(p(j)dx = --— \\/^(j)\^dx, 
2 dt J^3^ 

which relies on the fact that = on {xa = 0}. This proves our claim. □ 

4.2. estimate for the suitably linearized equations. We establish here an estimate for 
the solution {h,u,(j)) of the following linearized system: 

' dth + {ua + n)Vn + (n^ + n)V ■ ii + ii- V{na + n) + ndiv {ua + u) = rn, 



(4.9) 



^7 Tl ( Tl 

dtii + {ua + u)-Vu + u- Vua + T ( ) ) = V0 + r„, 

' Ua + n na \na + n' 



^ e^Act) = n + e "^"(/.(l + /i(0)) + r^. 



where r = {rn^fu^r^) is a given source term, and where h G {/io,/ii}, c/ (|4.6p - (|4.7p . The reason 
why we shall consider these two possibilities for h will become clear in view of paragraph 14.3.31 We 
add to the system the boundary conditions 

(4.10) '"31x3=0 = 0, (t)\x-j.=0 = 0. 
The crucial estimate is given by 

Proposition 2. Let {na,Ua-,4>a) the approximate solution constructed in Theorem \^ and some 
smooth (n,u,(j)) such that ^31x3=0 ~ '^'^^ 

(4.11) na + n>l/M, e-"^" (1 + /i(</>)) > 1/M, |n| + |n| + |(/>| < M, G [0, T], x G M^. 

Then, there exist C{M) and C{Ca,M) independent of £ (Ca only depends on the approximate 
solution) such that we have on [0,T] the estimate 



\\{h,u,(j),£V(p){t)\\l2(j^3^) < C{M)(j\{ho,^ J \\{e V^, 5tr<^, r„, r„) H^j^j^g 

+ C{Ca,M) (1 + ||Vi,a;(n,u,0)||^^^jj3_^ +e"^||n||ioo(K3_))||(n,n,(/),eV(/))||^2(^^ 



,) 
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Note that this is indeed a stability estimate for the Hnearized equation since when we take 
{u, n,(p) = that is when we hnearize exactly on the approximate solution, then we get from the 
above result the estimate 

^nj ^u) Ili2 

and hence from the Gronwall inequality, we obtain 

||(n,u,<^,£V0)(t)||i2 <e^'^*(||(no,tto)|li2+^ O^r^, r„, r„) H^^) , VtG [0,r] 

which is an type estimate for which the growth rate is uniform in £ G (0, 1]. 

Proof. In order to get this estimate, we shall use a linearized version of the total energy of the 
system. At first, let us collect a few useful estimates that we shall use for {na,Ua)- In the proof, we 
shall denote by Ca a number which may change from line to line but which is uniformly bounded 
for £ G (0, 1] and T G (0, Tq] where Tq is the interval of time on which the approximate solution is 
defined. Since the leading boundary layer term of vanishes, we have 

(4.12) sup \Ua\ + \Vx,tUa\ < Ca, Ca > 0. 

{0,T)xR\ 

For Ha, we have 

(4.13) sup \na\ + \Vxi,X2,tna\ < Ca, sup {dsflal < -Ca- 

To make more precise the last estimate, we observe that 

(4.14) dsna = ds (iVo {t,y, ^)) + 0(1) = ^dzNo (t,y, ^) + 0(1), 

where Z stands for the fast variable xz/e and hence, from the exponential decay of the boundary 
layer, we get that 

(4.15) sup 1 x353A^o(^,y,^3^ < \zd,N\t,y,z)\ < Ca 
te(o,T),a;eiR+ ^ V £ / te(o,r),j/eK2,2eiR+ 

Note also that ^a shares similar bounds: 



(4.16) sup \^a\ + |Va;i,a;2,t</'a| < C'o, SUp l^a^al < -Ca- 
(0,T)xK3_ (0,T)xR3_ £ 

with 

(4.17) 530a = ^3 ($0 (t, y, y ) ) + 0(1) = \dz^o [t, y, ^) + 0(1), 
and again, from the exponential decay of the boundary layer, 

(4.18) sup 1 X393$°(t,y,^) < sup \zd,N''{t,y,z)\ < Ca, 
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Let us now prove the energy estimate. First, multiplying the velocity equation by (n^ + n) ii, 
and performing standard manipulations, we obtain: 



(4.19) — / {na + n)-—= / {na + n)u-dtu+ / dt{na + n)—- 

(4.20) =Ii + l2 + h- I {ii-Vua) ■ ina + n)u+ (vu ■ {{na + n)u)) + dtiua + n 

Jr3 JjjS 



where 



h:=-ri ).((n,+n)n), h := i V<A • (n„ + n)n, 



Ua + n Ua \na + n 



[{ua + u) ■ Vn] • (ua + n)u. 

The last three terms at the r.h.s. of (j4.20p can be easily estimated by using ()4.12p . ()4.13p and 
(BUI: 



, {ru-{{na + n)u)) < C{Ca, M)\\ru{t)\\L2^^3j\u{t)\\L2^^s^), 
(4.21) £^ dtina + < C{Ca, M) (1 + \\dth\\L^) ||n(t)||i.(j,^), 

{U ■ VUa) ■ {Ua + n)u < C{Ca, M)\\u{t)\\l, y 



+ 



Let us turn to the treatment of /i. Integrating by parts, we first have: 

Vn \ \ \ f 1 // \ \ f V(na + n) 

■{{na + n)u)= / ■ — div ((na + n)u) - / • — -K-{na + n) -nu 



and hence, we obtain that 

(4.22) h = T / div ((na + n)^) + T • n u := /i + /i . 

yM3_ na + n J^3_ \ Ha Ua + n J 

To estimate we observe that 

Vn^ V(na + n) Vn^ n — UaVn 

Ua ria + n UaiUa + u) 

Consequently, we can use ()4.12p . ()4.13p and (j4.1ip to get that 

(4.23) < C{Ca, M){\\Vn\\L^ + e-^\\n\\L^) MlHrD W^^L^iK)' 
To estimate ll, we observe that we can write the first line of ()4.9p under the form 

(4.24) dth + div ((n^ + u)h) + div ((n" + n)n) = r„. 
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By using this equation to express div + n)uj, we obtain 

{dth + div {{ua + u)h) - r„) 



n 



3 n^j + n 



-a 



n 



1 



2(na + n) 2 



+ t: I n2(at + K + n)-V 



1 



'71°- + 



(4.25) 



1 f div (n" + u) 



+ / r„ 



n 



+ ria + n 



In the above expression, for the two last terms, we use (j4.1ip to get 



1 f 'h? div {ua + u) 



Ua + n 



Next, we observe that 



+ / Tr, 



n 



< CiCa,M){\\r„{t)\\L2 \\h{t)\\L2+a+\Nu\\L^) WmWl^ 



n 



dtiua + n) + (Ua + u) ■ Viua + n)] . 



From the equation satisfied by Ua and (j4.12p . we get that 

\dtna+Ua ■ VUal < Ca- 

Also, by using that {u)^ vanishes on the boundary, we have |(m)3| < 2;3||Vn||L;^ and hence by using 
(|4.15p . we also have 

(4.26) \\u-Vna\\L2<Ca\\Vu\\L^. 
Consequently, we obtain that 



+ n' 



< C{Ca,M){l + \\Vt,ML^ + \Nu\\L^)\\n\\l2. 



We have thus proven that 
(4.27) 

^ d r r*n2 

^ dt Jj^i 2{na + n) 



< C{Ca,M)(^{l + \\Vt,xn\\L°- + \Nu\\L-o)\\h\\h + {e~^\\n\\L-o + ||Vn||Loo) ||n||^2 Wuh^) + IKWh. 
As regards I2, since 



4>div {[ua + n)n), 



we use once again (j4.24p to write it as 



(4.28) 



(t)dth + j </)div [{ua + u)n) - j rn<j). 
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By differentiating with respect to time the Poisson equation in ()4.9p . we can express dtn in terms 
of (j), and substitute into the first term at the r.h.s of (j4.28p : 



4>dtrd>- 



For the last two terms, one has by using (j4.16p and (j4.1ip the straightforward estimate: 
(4.29) 

dt{e~^-{l + h{<t)))m'' - I <Pdtr^ 



< C{Ca,M){l + \\M\\l^)MI,^^^, + \\dtr^\\mRl)MLHRl)- 



As regards the second term at the r.h.s. of (j4.28p . we integrate by parts and use the Poisson 
equation to express n in terms of (j). We get 



(pdiv {n{ua + u)) = - I • {n{ua + u)) 

= - I V<j)-{s^A<j){ua + u)) + [ V<})-{e-^'^{l + h{cj)))<j)iua + u)) 

+ / V(j)- {r^{ua + u)) = Ji + J2 + J3. 



One has 



(4.30) J2 = J /" V|</.p • (e-'^'^il + h{cP)){ua + u)) = -\ ! div [e-^-{l + hm{'^a + u)). 

Once again, one has to be careful due to the boundary layer part of 4>a- Nevertheless, proceeding 
as for estimate (|4.26p . we obtain that [ua + u) ■ V0a is uniformly bounded in e and hence, we find: 

J2 < C{Ca,M){l + \\Vcl)\\L^ + \\Vu\\Lo^)U\\l2. 

Straightforwardly, thanks to (I4.12P one has also 

^3 < ^(Ca + ||w||l-) |k0||i2(R3_) ||eV(/>||i2(K3_). 

Finally, we compute 
(4.31) 

Ji=e^[ [{V^-V){ua + u))-V^ + e'' f ( ) • (n, + ^z) 

< [ \V^f\V{ua + u)\ - [ \^^diY (Ua + u) < C{Ca,\NxU\\L^)e^ [ |V</>|2 
where the last bound comes again from (I4.12p . 
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Combining the previous inequalities, we get the following bound: 
(4.32) 

^^ + \tAI e'|V0p + e'^'"(l + /i(</.)),^2 
2dt\ J^-i^ 

<C{Ca,M)[l + \\Vt,x(t>\\L^ + l|V,n||i«=)(||</>||2, +e2||v</,||2^) +e-2||^^||2^ + ||air^||2^ + ||^„||2^. 
Finally, to estimate Is defined after (j4.20p . we write: 

h = -\ {ua + n){ua + u) ■ V —- = / dvv {{na + n){ua + u))—-. 

2 Jlg3_ I 

Relying once again on ()4.12p and proceeding like for (j4.26p . we infer that: 

(4.33) Is < C{Ca,M){l + \\V{n,u)\\L^)\\u\\l2^^^^y 
Eventually, combining (OOll with (li:2TI) - (li:271l - (l02]l - (|i33]) . we obtain 

dt J^i \ 2 2{na + n) 2 2 / 

(4.34) + ^ \ . 

< C{Ca,M){l + \\Vt,Au,n,cl))\\L^ +e"^n\\L^){\\u\\l, + ||n||i2 + Vc^-lH^ + H^Hi^) 

We end the proof by integrating in time and by using (j4.1ip . □ 

4.3. Nonlinear stability. We shall now work on the nonlinear system (j4.ip in order to get The- 
orem [3l Thanks to the well-posedness in H"'' for m > 3 of the system ()4.ip . we can define 

= sup{r e [0,ro], Vt G [0,r], ||(n,u,0,eV0)||^^„(R3) < e^ and gSD is verified} 

where r is chosen such that 

(4.35) 5/2<r <K 
and the if™ norm is defined by 

\\f\\HriRi)= E ^'"'ll^"i^"2^S/IL2(R3^). 

|Q![<m 

We shall also use the norms: 

(4.36) ||/(t)||^^^_^(j,3^)= Yl w^o^^r^r^rfmL^iRi) 

\a\<m 

where the vector fields Zi are defined by 

(4.37) Zo = edt, Zi = edi, i = l, 2, Z3 = ^3 

1 + a^s 

and 

ll/WllHrK)= E £'"'ii5r9.",^5.",^5,"33/(t)ii^2(i,^). 

For the sake of brevity, we will also use the following notation: 

Z" = 1 ^2"' ^3' ' for a = (ao , ai , 02 , as ) • 

Finally, we set 

Qm{t) = ||(n,'u,0,eV0)(t)||^m_^ + ||w(t)||^m-i 
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where we have set u = ecurlu i.e: 

/ d2U3 - d3U2\ 

(4.38) uj = e I d^ui - dms . 

\diU2 - d2Ui J 

We shall first prove that Qm is the important quantity to control for the continuation of the 
solution and then we shall estimate Qm{t) by using Proposition [2j 

4.3.1. Estimate of the T-L^ norm. Let us recall a classical estimate for products in dimension 3: 

(4-39) lk^llL2(]Rp < ||w||^n(lg;3_)||?;||^S2(]R3 ), 

with si + S2 = 3/2, si / 0, S2 / 0. More generally, 

(4-40) ||ni . . .uk\\l'^(^^3^) < ||^ii||//n(R3_) . . . \\uk\\h'k{m?^) 

for si + • • • + = i(K - 1), si / 0, . . . , Sp 0. 

As an application, we state 
Lemma 1. For any u, v and e G (0, 1], we have uniformly in e 

(4.41) \\{edt)Huv)\\Hi^Ri) < \H\L^\\{£dt)''v\\Hi^^3^) + \\v\\L^\\{edt)''u\\Hi^^3^) 

with the notation {edt)^f = {f,edtf,--- ,{£Ot)^f) if k > 0, if k < 0. Moreover for any smooth 
function F with F(0) = 0, we have 

(4.42) ||(e5i)'=F(n)||^,(K3) < C (|| (e9i)^||^,(K3^) + \\{edt)'-\\\Hl+\B.l)) ^ ^ > 1' 
with 

C = C[||u||i^(R3_),e-t||(eat)'=n||^,(K3) + e"^IKe9t)''"^'"llHi+i(R^)]- 

Proof. Let us first emphasize that a simple rescaling {t' = t/e, x' = x/e) allows to restrict to the 
case e = 1. 

To prove ()4.41|) . it suffices to use the Leibnitz formula. When all the derivatives are on u or 
on V, we estimate the other term in L°°. For the remaining terms which are under the form 
d^'d^^u {dtf^d'^-^v with /ci + |ai| < A: + / - 1 and A:2 + |a2| < + / - 1, we use ()T39]) to write 

with si + S2 = 3/2. Therefore, by taking si and S2 smaller than 1, we obtain ()4.4ip . 

The proof of (|4.42p (still for £ = I) relies on the Faa Di Bruno formula for the quantity d^d^F{u), 
\a\ < I, and on the product rule (|4.40p . For brevity, we just treat the case 1 = 1. 

If |q| = 0, dfd^F{u) = dfF{u) can be decomposed thanks to the Faa Di Bruno formula as a 
sum of terms of the type 



F^P\u) f[ d^'u, with Y.ki<k. 



i=l 



Either there is zero factor in the product (that is /c = 0), and we use the bound 

II^(^^)IIl2{R3) = ll^(^^) -^(0)||l2(k3_) < C[||n||Loo] ||n||i2(R3_) 

or there is one factor in the product, which leads to 

l|i^'H5M|i2(M3) < C 
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\(^t ^IIl2(R3 ) 



or there are several factors (which means ki < k for all i), and we have by inequality (j4.40p (with 
si, . . . , Sp-i close to 3/2 and Sp = 1) 

p 



\\F(PHn)lldtuh.^ 



L) 



< c 



\u\\L'^{Rl)A\{dt) ^ll//3/2(]g3 ) 



i=l 



If |a| = 1, for instance a = (1, 0, 0), d^d^F{u) can be decomposed as a sum of terms of the form 



F(P+i)(n) Yl dS (dfd^.u) , with + ^" - 



i=l 



Either k" = k, and we use the bound 

or k" < k — 1. In this latter case, we apply again (j4.40p : if one of the k'^'s is k, we get 
or all /cj's are less than k — 1, so that 



\\F^p-''\u)lldS (afa.,«) 11^2(^3) < c \\uU^,\\{dt)'-\\\HSf2 



i=l 



The last inequality comes from (|4.4U|) with si = ■ ■ ■ = Sp close to 3/2 and Sp+i 
dt" dxiu). Combining above inequalities is enough to obtain ()4.42p . 



1 (for the term 



□ 



We shall first prove that by using the equation, we can estimate the T-L^ norm of the solution of 
m on [0,r^) : 

Proposition 3. For m> 3, we have for every e £ (0, 1], for t G [0, T^) 

(4.43) \\{n,u,<P,eV<P){t)\\L^^^3^^ < Ce'-^/^, ||V(n,n,</),eV</))(t) 1^00(^3) < Ce^-^/^ 
for some C > independent of e and 

(4.44) ||(n,^z,0,eV0)(t)||^^(K3^) < C[a, M] (e^+i + e^) 

where C stands for a continuous non- decreasing function with respect to all its arguments which 
does not depend on e. 

Proof. The first set of estimates can be obtained by using the Sobolev inequality in dimension 3 
and the definition of . 

To prove (j4.44p . we proceed by induction on the number k of time derivatives. For A; = 1, we 
proceed as follows. At first, by using the evolution equations on n and u in system (|4.ip . we can 
compute edt{n,u). We claim that, for t S [0, T^) 



(4.45) ||e9t(n,u)(t)||j^,n-i(jj3^) 



< C 



Ca,M,\\{u,n) 



l°°(ir3)> lkV(u,n)||^cx>(K3),e ^ "')ll/f™(R3 



Indeed, the expression of edt{n,u) involves four kinds of terms (besides the "easy ones" eV0 and 



i) terms that are linear in e(Vn, Vn), with coefficients depending on Ua,na- They are bounded 
by Ca ||eV(n, (i;3 ■) ) where Ca depends on the norms of (eV^;)" (lia, n^), with |a| < m—1. 
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ii) terms that are linear in {n,u), with coefficients depending on najS'Vna, Vn^- They are 
bounded by C[Ca,M] ||(n,n)||j^m-i(]jjp. 

iii) quadratic terms, involving products of n or n with eVn or eVn. By ()4.4ip . Lemma [H we 
can bound them by 

3 

C ||(n,n)||Loo + ||eV(n,n)||i^(jj3_) +e"5||(u,n)||j:^™-i(jg3) 

iv) a fully nonlinear term, coming from the pressure in the Euler equation for the ions. It reads 

ln(na + n) ln(na + n) 



eVF{na,n), with F{n a, n) := T 



ria + n 



n„ 



In particular, F(na, 0) = 0. The nonlinear term can then be evaluated in using (j4.42p (with 

k = 0,1 = m). It is bounded by 



C 



n 



Combining the previous bounds yields (j4.45p . Therefore, by using ()4.43p and the definition of 
T^, we obtain 

(4.46) ||eOi(n,n)(i)||^^™-i(R3) < C [c„ M, e-^-^/s] (g^+i + 
Moreover, by applying edt to the Poisson equation in (14. ip . we get that 

{e^A - e-^^"+^^)edt^ = edtn + edt{e-'t'-) (e""^ - l) + e'^+^edtR^, edt<P/.,,=o = 
and by applying next {edx)", 

(4.47) (e^A - e-('^"+'^)) {edx)''edt(l) 

= {edxredtn + {edxrGi{-,(t>) + [{edxr;G2{-,4>)]edt4> + e''+\edxredtR^, 

where Gi(-, (/)) = edt{e~'^'') [e~'^ — l) , G2 (•,</)) = e~^'^'^~^'^\ From there, one can perform standard 
energy estimates, recursively on |a|. The nonlinearities are handled thanks to (I4.42p . This leads to 



(4.48) ||eat(0,eV0)||^^^-i(j,3) <C 



Ca,M, ||(/)||^oo(iR3_),e ^^"^WwH^iR^) 

(Ik5t^ll//--1(R3) + ||(/'||j:^™-l(R3) 

Consequently, by combining the last estimate and (I4.46P and by using that r > s, we get that 

||e5t(n,n,</.,eV0)||^™-i(jj3) < C[Ca, M]{e''+^ + e^). 
Now let us assume that we have proven that 

(4.49) ||(eai)'=(n,n,(/),eV</.)||^™-.(j,3) < M] (e^+i + e^) . 

By applying e (edt)^ to the evolution equations for n and u in (14. ip . we obtain an expression for 
{edt)^+\n,u). This expression invoves k e-derivatives with respect to time of linear, quadratic 
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and fully nonlinear terms. These terms are evaluated again thanks to Lemma [TJ We get that for 
k + l<m, k>l 



\\{edt)''^Hn,u)m^^... 



-) 



< C 



Ca,M, ||(n,n)||i«,(iR3), ||eV(n,n)||i^(jR3),e ' ^ \\{edt)'' ^ {u,n)\\j^^-k+j , 



j=0 



+ e 



K+l 



and hence from the induction assumption, we obtain 

(4.50) \\iedt)''-Hn,u)it)\\^.^-,-,^^,^^ < C[a, M] (e^+i + e''). 

Finally, we can use the Poisson equation to control {edt)^^^{e'S/(j), (j)). More precisely, applying 
{edt)^~^^ to the Poisson equation, we get 







where C''{(f)) is a commutator. Thanks to Lemma [Hand the induction assumption, we have for this 
commutator the estimate 



< c 



j=0 



Y.\\{edt)'"^<P\\^.^ 

j=0 



m — k / 



Consequently, the standard a priori estimates for the elliptic equation yield 

and hence we get by combining (|4.50p and the last estimate that (j4.49p is verified for k changed in 
k + l. 

□ 

By using similar arguments as above, we shall also get: 
Lemma 2. We have for e G (0, 1], the estimate 



(n,n, (A, eV (A)(0) ||^™(r3^) < C[Ca,Ml 



.K 



Proof. The proof follows the same lines as the previous Proposition. The main difference is that 
because of the choice of the initial data, the estimate (j4.46p is replaced by 



The end of the induction can be performed in the same way. 



□ 
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4.3.2. Normal derivatives estimates. The aim of this subsection is to prove that we can replace one 
normal derivative by tangential derivatives thanks to the equation and hence that ||e ^3 (n, ^3) (t) ||^m- 

can be estimated in terms of Qm- We shall again need some product estimates: 
Lemma 3. For any u, v and e G (0, 1], we have uniformly in e 

(4.51) ||(ea3)'=^°MllL2(M3^) < lkl|L-||(e93)'^||^H^(K3) + ||HlL-||^ 

Moreover for any smooth function F with F{0) = 0, we have 

(4.52) ||(ea3)'=Z"F(n)|L.(^3) < Ci\\{eds)'^u\\^,^^j^^^^,^^ 

1<I/3|<2 



where 



Ci = Ci 



Ci = C 



_3/2|| II 

-)'^ ll^ll^fe+l«l(R3 ) 



are two continuous, non- decreasing functions independent of e. Moreover C2[',0] =0. 

Note that the last statement of the lemma implies that C2 is small when its second argument is 
small. This will be used in the sequel to absorb some error terms in the estimates. 

Proof. To prove (j4.5ip . it suffices to use again the Leibnitz formula. When all the derivatives are 
on u or on u, we estimate the other term in L°°. For the remaining terms which are under the form 
{ed3 f^Z'^^u{ed3)''^Z"^v with fci + |ai| < /c + |a| - 1 and k2 + |a2| <k + \a\-l, we use KMi to 
write 

with si + S2 = 3/2 and we get the result by taking si and S2 smaller than 1. 

The proof of (|4.52p is very similar to the one of (j4.42p : broadly, (ed^) substitutes to eV, whereas 
Z substitutes to edt. The full gradient terms (eV)^, |/3| = 1,2, are still connected to the use of 
(|4.40p . with Si £ [1,3/2]. We leave the details to the reader. □ 

With this technical lemma in hand, we now prove the 
Proposition 4. There exists eq such that for every e G (0,eo]) t £ [0,^^); < /c < m, 
||(e53)'=(n,n,</.,eV</.)(t)||^™_fe(ig3^) < C[C„ M] (e^'+i + Q^(0). 

Note that we can reformulate the above Proposition into 
(4.53) ||(n,7x,0,eV0)(t)||^^(K3^) <C[Ca,M](e^+i + Q„(t)), VtG [0,T^). 

Proof. By using the definition of the vorticity, we already have: 

||ea3(ni,n2)|Lm-i < ||t^|Lm-i + < Qm{t). 

Next, we notice that 
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By using the Poisson equation to express e^d^'^cf), together with Lemma[3]to control the nonhnearity, 
we get that 

,£-3/2 I 



\\ed,{eVm\\ 



njm — l ^ C 



< C 



Finally, by using the equations on n and U3 in ()4.ip . we get that 
(4.54) 

, ( dsns \ _ 

/ dtn + {ua + n)i,2 • Vi,2n + (1 - x){ua + u)^d^n + u ■ Vua + ndiv ti^ + (n + na)Vi,2 • {ui,U2) - Rn 

I dtUz + (Ua +'u)l,2 • 

where An is defined by 



Note that thanks to ()4.6p . we get that A„ is invertible, its inverse being given by: 



A. 



-1 



1 



T^(n" + n)-i -xK + n)3 

-X(tia + u)z Ua + n 

Hence, we can use this system to estimate ||e93(n,ii3)(t)||^m-i^jg3 y Note also that the field e(l — 

x)c^3 is equivalent to Z^. One can again decompose the terms into 

i) linear terms, whose coefficients depend on Ua,'^Ua,na,£Vna and involve {u,n) or Zi{u,n). 
They are bounded by C[Ca, M]\\{n,u)\\u^^^^. 

ii) nonlinear terms: they all involve products of the type F{x,n,u) Zi(n,u), where the function 
F satisfies F{x, 0, 0) = 0. They can be estimated with Lemma [3] (A; = 0, |a| = m — 1). 

We obtain 



(4.55) ||ea3(n,U3)(t)||^^^-i(jj^) < Ci (\\{n,u)\\n^^^^ + W^^^uToTeH^D + ^^^^ 

+ C2(\\in,u)\\nz,, + Il(e53)(n,w)||^™- 



with 

and with 



Ci = Ci 



Ca,M,\\{n,u) 



\Z{n, u) 



Co — Co 



Ca,M,\\{n,u) 



\Z{n, u) 



-3/2 1 



that vanishes with its last argument. We used here the notation Zf = {Zif)Q<i<3. 

Now, from the L°° estimates and (I4.44|) in Proposition [3] and from the definition of T^, it follows 
that 



|e53(n,U3)(t)||^m-i(R3^) < ^^^1 (||(n, u)||^„ ^(^3^) + ||eV(/.||^™-i(^3 ) + e' 



where C2 vanishes with its argument. As r > 3/2, we can absorb the last term at the r.h.s in the 
l.h.s for e small enough: this yields 

(4.56) ||ea3(n,n3)(t)||^^^-i(^3)< (||(n,n)||w,n_^ + ||eV</)y^ + e^+i) < e^+i + Q„(t). 



I 'Tjm— 1 
' lico 

We can then prove Proposition [H by induction on k. 
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The estimate for /c = 1 is thus already proven. Let us assume that it is proven for j < k. At 
first, from the definition of the vorticity, we find that 

\\{ed3)''^\ui,U2)\L^-(k+i) < \\{ed3)''uj\\ m-(k+i) + \\{ed3)^u\L,n-k 

't-co,£ 't'Co,e ^co,e 

and the right hand side is already estimated thanks to the definition of Qm and the induction 
assumption. In a similar way, by using the Poisson equation, we get that 

k 

and we conclude again by the induction assumption. Finally, to estimate \\{ed3)^^'^{n, U3)\\ m-(k+i) , 
we use again ()4.54p . Lemma [3] and the induction assumption. 

□ 

4.3.3. Main energy estimate. We shall finally estimate Qm by carefully using the system (j4.ip and 
the linear stability estimate of Proposition [2j 

Let us before state a useful commutator estimate: 

Lemma 4. For i = 1, 2, 3, we have the estimate: 

(4.57) ||[^",/]a,<7||i2(K3) < C[||(Vt,,./,a,5)||ioo(R3),e-t||/||^i.i(^3^)] 

(Il/ll^hl(]g3_) + lbll^l"l(K3))- 

Proof. To get this lemma, it suffices to combine the Leibnitz formula and (j4.39p as in the proof of 
Lemma [Hand Lemma O The only term that we handle in a different way is when all the derivatives 
but one are on dig, in this case, we write 

□ 

Our main energy estimate for ()4.ip is the following: 
Proposition 5. There exists eo > such that for every e £ (0,eo)) have the estimate 

QUtf <C{Ca,M){e^^ + Te'^ + /*Q^), yt e [0,T^). 

Jo 

Proof. We shall first estimate \\{n,u,(j),eV(j))\\-^m. (m.\)- To this end, we apply the operator Z" to 
the system (j4.ip . We obtain: 

' dtZ^n + {ua + u) ■ VZ^'n + {ria + n)div Z^'u + Z"'u ■ V{na + n) + Z^'n div {u + Ua) 

~ Cn ~\~ S Z 

--^l 

Ua + n Ua \na + nJJ 

= VZ"(j) + Cu+e^Z''Ru, 
e'^AZ'^^ = Z'^n + e-'^»^"0(l + h)+C^ + e^+^Z'^R^. 

One has h = ho for |q| = and h = hi for \a\ > 1. The functions Cn,Cu and are remainders 
due to commutators. One can observe that this corresponds to the "abstract" system that we have 
studied in Proposition [2j We shall now estimate the remainders in order to be able to apply the 
stability estimate ()4.34p . 
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We first claim that 

(4.59) \\{Cn.Cu)\\L^ < C[Ca,M,e'-'-^]\\{n,u,^,eVmHT(^l) ^ C[Ca. M]{e''+^ + Q^{t)) . 

Tliis estimate is in the same spirit as the previous ones. One must distinguish between linear terms 
and nonlinear terms, the latter being controlled thanks to the product estimates of the previous 
lemma. For brevity, we only point out two typical terms: 

• a typical linear term in Cn is [Vn^; Z"'] ■ u. It is a priori singular (because of the boundary 
layer), but noticing that 

\\[Z,Vna]f\\L^<Ca\\f\\L^ 

gives the bound \\\SIna,Z'^] ■ n||^2(K3) < Ca ||u||^m(i.3_) < C[Ca,M] + Qmit))- 

• a typical nonlinear term in Cu is the commutator 

r [F(n„,n)V;Z"]n = r [F(na, n); Z"]Vn + rF(n„, n)[Z"; V]n 
with F{na,n) := ^ - Clearly 

||F(na,n)[2:'*; V]n||i2(iR3_) < CM\\n\\^\c,\^^2^^ 
As regards the first term at the r.h.s., we use the commutator estimate of Lemma [U 



||[F(n„,n);2:"]Vn||^2(R3 ) <C || Vt,xi^(na, n), V^.n||^cx,m3), e ^/^||F(na, 



n 



)^'^ II vt,x-tH"'a' '^i'"-|lL°°(R5.))^ IK "^ll-H; 

(^||F(na,n)||^™.(lB,^) + ||n||^m(iR3_)y . 

The product estimate (I4.42P implies that 

ll-^(«a,n)||^^(R3_) < C Ca, M, ||n||^cx,(R3_),e"^/^||n||^^(iB,: 
so that eventually 

\\T [F(n„n)V;Z"]n||i2(K3^) < C7[C„ M, e'-5/2] ||n||^„(K3 ) < C[Ca,M] {e^+^ + Q^{t)) . 
The treatment of the other terms is left to the reader. 
For the commutator C^, we get that 

||C^||l2 < C[Ca,e''"^/^]e||0||^m-i(]g3) + e (||eV(^||^m^^(K3) + We"^ HnT'^^l)) ' 

The first term at the r.h.s corresponds to the commutator with the semilinear term in the Poisson 
equation. The second one bounds the commutator [e^A;^"]^. We use the fact that 

[^3, ds] = -83 ( — ^ ) £^3. 



Consequently, by using again Proposition [3] and Proposition U we get that 

(4.60) \\C4l^ < CiCa, M, e^-3/2) (e^^+i + Q^) . 
By using similar arguments, we also get that 

WdtC^L^ < C(C7„, e^-3/2)||^||^„(jj3^) + ||e2v20||^,„^^(jj3) + \\eV ^^n. ^^s^y 

To estimate the second term, we use the elliptic regularity for the Poisson equation in (14.58p . This 
yields 

and hence, we find that 

(4.61) WdtC^L^ < C7[C„M,e'-3/2](^^+i + g^). 
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.K+1 



Note that thanks to Proposition [3l we have that on [0, T^), 



lL°°(R'^) 

and hence that the assumption 14.11] in Proposition [5] is matched on [0,T^) for e sufficiently smah. 
Consequently, we can use Proposition [5] and Lemma [5] for the system (j4.58p to get 



2 

L2(R3 ) 



|2 



||(n,u,(/>,eV(^)(t)||^„ ^^jj3) < C[Af](^e^'^ + / ||(e ^r^,dtr^,rn,ry) 

+ C[Ca,M] / (l + \\Vt,x{n,u,(j))\\^^^^s^ +e"^||n||^oo(R3_))||(n,n,(/>,eV(?;»)) 
Jo ~^ 

where we have set 

rn = Cn + e'^Z'^Rn, ru = Cu + e'^Z'^Ru, r^ = C^ + e^+^2"i?<^. 
Consequently, by using the estimates ()3.7p . ()4.59p . ()4.60p . (j4.6ip and Proposition [3l we get that on 

(4.62) ||(n,7x,0,eV0)(t)||5^„^(^3) <C(C„M,e'-5/2)(e2if^^2i^^^ j Q^itf^ 

It remains to estimate ||a;||.^m-ij.jg3 ^ i.e. the second term in the definition of Qm- By applying 
the operator eVx to the second equation in the system (|4.ip . we find that 

(4.63) dtuj + {ua + u)-Vuj = er{ua,Vua) (Vn, Vn) - eu ■ + R^ 

where V{ua, Vn^) is a polynomial of degree less than two whose coefficients depend only on Ua and 
Vti(j and thus are uniformly bounded. Note that we have used that 

/ \/ fi \/ fi / fi w 

Vx — -[ — — = V X V(log(na + n) -logn«) =0. 



^ria + n Ua \na + nj 
By applying the operator (e3)° = {edt^x)" to ()4.63p . we find that 
(4.64) dt{ed)''uj + {ua + u) ■ V{ed)''uj = C 

where by using Lemma (TJ we have 

||Ck2(M3) < C[C„M,e'-5/2](||^,||^^(^^) + ||[(ea)";^x. V]^||i2(K3^) +e^). 

To estimate the commutator, we use the following variant of Lemma H) we first write since 
||cij||T,m-i < that 

II \ \ He ^11 II n.j- 

||[(ea)";n- V]a;||i2(K3) < ||Vj,:,n||ioo(R3_) ||'u||^^(r3_) + || (eS)""?! • Va;||i2(R3 ) + e"^ ||u||^„(jj3 ) 
and to estimate the second term, we use also (j4.39p . to get since |a|<m — l,m>3 that 

||(e9)"u • Vwll < e"^||u||wm(K3)||Vtj||^i(]K3_) < £~^\\uf^^^^A^y 
We thus get the estimate 

l|Ca;||L2(]R3_) < C[Ca,M,e''"^/^]||'u||^„(]j3_). 

Consequently, from a standard L? type energy estimate on the equation ()4.64p (let us recall that 
{ua + u)^ vanishes on the boundary), we get that 



-IImI|2 



<C7[Ca,M](e2^ + ^Vll?^^(R3)) 
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and hence by using Lemma [2] and Proposition HI we infer that for t E [0, T^), 

(4.65) ||l^||5^™_1(^3^^ < C[Ca,M]{£^'' + e^''t + I^Ql). 

To end the proof of Proposition El it suffices to combine the estimates (|4.65|) and (|4.62|) . 

□ 

4.3.4. Proof of Theoreml^i We can now easily conclude the proof of Theorem [3l By using Propo- 
sition [5] and the Gronwall inequality, we get that 

Qm{t) < C[Ca, M]e^e^[^-^-^l*, Vt G [0, T'). 

This yields thanks to (j4.53p . 

\\in,u,cP,eVcP)it)\\Hrn^^3^) < C[C„M](e^[^-^l* + £)£^, Vt G [0,r^). 

Note that we also have by Sobolev embedding 

||(n,0)||ioo(M3^) <e^-3/2(e^[^-*-^l* + e), VtG [0,r^) 

and 

||X(X3)^3|L.(I,3^) < C[J].^-^/2(e^[^-^l* + e), € [0,T^). 
In view of the definition of T^, we obtain that for e sufficiently small > Tq and that 

\\{n,u,^,eV^){t)\\Hr.(^^z^) < C[Ca,M](e^[^-^^l* + e)e^, Vt € [0,ro]. 
This ends the proof of Theorem [3l 

5. Further remarks 

5.1. Generalization to general smooth domains. In this paragraph, we briefly explain how 
the work which has been done for the half-space Mi]_ can be adapted to handle the case of a smooth 
open set Q C M']_ whose boundary dQ is an orientable compact manifold. 

i) Derivation of boundary layers. 

The principle is to apply the tubular neighborhood theorem in M^; then, dQ being orientable, 
there exist a smooth function (p : dQ — )• M (normalized to have \Vip\ = 1) and a neighborhood 
U := {x G il., if < r]} of in il. such that for any x £ U, there is one and only one (x, y) G 9$7xM+*, 
such that 

X = X + y Ux- 

where is the inward-pointing normal at x. 

Then we can look for approximate solutions of the form: 

K 

«pp,ulpp,(Plpp) = '^e' {n\t,x),u\t,x),(l)\t,x)) 

(5.1) 



+ E (iV* (t, i, , {t, X, ^) , <^>' (t, X, ^)) . 

i=o see 

and the construction of Section [3] remains the same, 
ii) Stability of the boundary layers. 

The stability estimate that is based only on integration by parts is still true. For higher 
order estimates, one can modify the definition of the conormal Sobolev spaces in the usual way by 
considering a finite set of generators of vector fields tangent to the boundary. 
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5.2. Generalization to non homogeneous conditions for u. Let us now explain how we could 
generalize our work to the following boundary condition on u: 

(5.2) ^31x3=0 = Ub. 

We shall impose a subsonic outflow condition —Vt^ < ut, < 0. One can readily check that the 
boundary {xs = 0} is noncharacteristic in this setting and that this boundary condition is maximal 
dissipative. From the physical point of view, this boundary condition means that some plasma is 
somehow absorbed at the boundary. 

Let us once again explain what are the main changes in each principal part of the proof. Details 
are left to the reader. 

i) Derivation of boundary layers. 

Considering the derivation of the approximation, computations get more tedious, but the results 
remain roughly unchanged. In particular, we can check that there is still no boundary layer term 
for the velocity at leading order, which means with the same notations as before, that = 0. 
The principle to show this fact is to obtain closed equations on Uy and t/g which have the trivial 
solution 0, as it was done for Uy in Section [3l 

ii) Stability of the boundary layers. 

The main part which is modified due to the new boundary condition is the proof of Proposition 
[21 in which new boundary contributions appear since {ua + u)s\x3=o 7^ 0: we have to treat them 
carefully. Nevertheless, we note that since the boundary condition for the original problem is 
dissipative, the boundary terms when we integrate by parts in the transport terms have the "good 
sign". For example, one typical new harmless term comes for instance in ()4.25p . By integration by 
parts for the term 

r V(n2) . (m"^ + u) 
Jul na + n 

the new term that appears is — J^2 ^^"^"^^'^"'^^^ 1x3=0 and one can observe that {u"" + u)^\x3=o < 0, 
so that this term has indeed the good sign. 

One more significant term comes from the computation of Ji in (I4.3ip . One can check that a 
first integration by part yields the term B := f^2 |V(/>p(na + 3 1^3=0 and a second integration by 
parts yields the term —^B, so that only the sum of these two terms, ^B, has the good sign. 

Furthermore, some terms also have to be estimated more carefully. In ()4.29p and (j4.30p . the two 
terms 

-If dtie-'^^il + him<P\'-l [ |0|'div(e-^'"(l + M0))K + n)) 

can not be estimated separately as before. They can be treated using the equation satisfied by 
e~'^'' (that is the transport equation satisfied by Ua, up to some error term). 

The end of the proof of Theorem [3] is unchanged, except for the last part which gets actually 
simpler. Indeed, the boundary {xs = 0} is now noncharacteristic and therefore we can estimate all 
normal derivatives using directly the equations on {n,u), without having to introduce the vorticity. 

When < Uf, < Vt^ the boundary {x^ = 0} is still noncharacteristic (this corresponds to some 
inflow condition, which means that some plasma is injected at the boundary) but we still need to 
add two (in dimension 3) boundary conditions in order to get a well-posed problem. There are 
many possibilities and some of them do not yield a dissipative problem. In this case, our approach 
to the stability estimate of Proposition [2] would break down. We leave this case for future work, 
as well as the case of supersonic outflow condition, relevant to plasma sheath layer problems. 
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5.3. The relative entropy method. In this paragraph, we focus on the one-dimensional case. 
First, we recall that Cordier and Peng constructed global entropy solutions to the same isothermal 
Euler-Poisson system as the one studied in this paper, posed in the whole line M (see [8]). Although 
it has never been done, it seems reasonable to believe in the existence of some global entropy solution 
in the half line M"*", at least for the simplest boundary condition 0;, = (and non-penetration, that 

is U3|a:3=o = 0). 

Assuming the existence of such solutions with weak regularity, one could study the quasineutral 
limit, using the so-called relative entropy method which was introduced by Brenier in [3]. This 
strategy was previously briefly evoked in It consists in considering the following modulated 
nonlinear energy: 



(5.3) 



'We(i) =^ j n\u - Uapp\^dx + j "-(log {n/uapp) -1 + napp/n)dx 
+ j (e-'^log (e-'^/n^pp) - + napp)dx + ^— j \d^(l)\'^dx, 



Then, by an explicit computation (similar to the one given in [21 j ) of the derivative in time of 
He-, the principle is to show that this is a Lyapunov function, and thus that it vanishes as e — )• 0, 
if 'He(O) — )-e->.o 0. This would roughly proves strong convergence in for {n^,Ue,(f)e) without 
performing high-orders estimates in order to obtain strong compactness, as in the proof of this 
paper. Notice furthermore that this seems anyway impossible to justify such high-order estimates 
for solutions with low regularity. 

The convergence would still be local in time, since for this method, one has to consider strong 
solutions (that is with at least Lipschitz regularity) to the limit systems. 

5.4. Two Stream instability. It seems reasonable to consider the case of several species of ions, 
say two for the clarity of exposure. This means that each species, described by (ni,ni) (resp. 
(n2,n2)) satisfies an isothermal Euler-Poisson system: 



(5.4) 




-I- div {uiUi) = 0, 

+ Ui ■ Vui + T' Vln(nj) = V0, 

and these systems are coupled through the quasineutral Poisson equation: 

(5.5) e^A(f) + e~'^ = ni + n2. 

One could expect to prove similar results to those proved in this paper. Nevertheless, the 
situation turns out to be very different in that framework (even without boundary). This is due 
to the so-called two stream instability. We refer to the work of Cordier, Grenier and Guo ^7j, who 
studied that mechanism and showed that any initial data such that 

is non linearly unstable, hence the name two stream instability (actually they do not explicitly 
treat the case of ()5.5p . but their study also holds in that case). Instability means here that there 
exists some ao > such that, for any 6 > , there exists an initial condition in a ball of radius 
6 (for a Sobolev norm with regularity as high as we want) and with center the two-stream initial 
condition and whose associated solution goes out of the ball of radius qq in the norm, after some 
time of order 0{\ log5|). In other words, there are perturbations of the two-stream data which are 
arbitrarily small in any Sobolev norm (as high as we want) and which are wildly amplified. 

One can then observe that it is possible to go from the system (j5.4p - ()5.5p with any e > to 
(I5.4p - (l5.5p with e = 1 thanks to the change of variables {t,x) i— )■ (fjf)- This means that the 
quasineutral limit can somehow be seen a long time behaviour limit for the Euler-Poisson system. 
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By extrapolating the two-stream instability results, one can infer that the formal limit for (j5.4p - (j5.5p 
is false in Sobolev regularity for such data. 

On the other hand, one can conjecture that the formal limit is true as soon as we avoid two-stream 
instabilities, that is to say when: 

■"1,0 = ■"2,0- 
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